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^ ABSTRACT 

We study the collective density modes that can affect neutron-star thermodynamics in the baryonic density range 
between nuclear saturation (po) and 3po- In this region, the expected constituents of neutron-star matter are mainly 
neutrons, protons and electrons [npe matter), under the constraint of (3 equilibrium. The elementary excitations of 
this npe medium are studied in the random-phase approximation framework. We emphasize the effect of Coulomb 
interaction, in particular the electron screening of the proton plasmon mode. For the treatment of the nuclear 
interaction, we compare two modern Skyrme forces and a microscopic approach. The importance of the nucleon 
effective mass is observed. 
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^ ■ I. INTRODUCTION 

>, 

\ The evolution of neutron stars is in general determined by the microscopic processes that can occur in their interior 
• and in their crust. In particular, the long-term cooling by neutrino emission is strongly affected by the detailed 
structure of the matter and the excitations it can sustain. The mean free path and emissivity of neutrinos are directly 
related to the correlations present in neutron star matter, which in the standard models is composed of neutrons, 
protons, electrons and muons [1]. The effects of collective excitations on neutrino emissivity have been extensively 
studied in the literature, often with controversial results [2-10]. The plasmon excitations seem to play a relevant 
role in quenching the contribution of the pair-breaking mode to neutrino emissivity [7] , since it obscures the possible 
effect of the Goldstone mode that should be present in a neutral superfluid due to gauge invariance. However, since 
^ ' neutrons, protons and electrons are coupled, it is not clear to what extent the presence of the plasmon mode can 
affect the whole set of collective excitations or the excitation spectrum at a qualitative level, even for nonsuperfluid 
matter. It has to be stressed that the plasmon is mainly an electron excitation, since the electron screening on 
proton-proton Coulomb interaction is quite effective, but the proton density oscillations can become relevant at not 
too low momentum. Furthermore, some of these excitations are damped and it is essential to estimate which actually 
are overdamped and therefore must be excluded. 

Collective modes in asymmetric nuclear matter have been studied previously, e.g. in Refs. [11-13]. In the astro- 
physical context, a study of the collective excitations in normal neutron star matter on the basis of the relativistic 
mean- field method has been presented in Ref. [14]. In this article, we present a study of the excitation spectrum in 
homogeneous neutron star matter under realistic conditions of proton fraction at different densities, paying specific 
attention to the role of electron screening and to the damping of the modes. Different nuclear models are used to 
see the dependence of the spectrum on the nuclear interaction: two modern Skyrme effective forces (SLy230a and 
NRAPR) are compared with effective forces deduced from microscopic Brueckner-Hartree-Fock calculations. In both 
types of approaches, the nucleons are considered nonrelativistic. We treat the nucleonic part of the response function 
with the quantal random-phase approximation (RPA), while the semiclassical Vlasov approximation is used for the 
electronic part (ultrarelativistic) . Since in this article we want to clarify the relevance of the nuclear and Coulomb 
forces, the infiuence of pairing and the problem of neutrino emissivity are left to future works. 

In Section II, we briefly remind the RPA formalism used to obtain the dispersion relations and strength functions. 
Section III is dedicated to the study of electron-proton coupling through the pure Coulomb interaction, to emphasize 
the suppression of the proton plasmon mode by electron screening. The realistic situation, involving Coulomb and 
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nuclear coupling among neutrons, protons and electrons, is analyzed in Section IV. Conclusions are exposed in 
Section V. 

II. RPA FORMALISM 

We will study the elementary excitations in homogeneous nuclear matter in the RPA, which amounts to considering 
correlations only at the level of the one-particle— one-hole interaction. 

Let us first consider a single fluid. The RPA equations for the polarization propagator Il{q,u)) is a function of the 
momentum q and energy w. 

U{q,w) = Uo{q,w) + t;res(5)no(g, a;)n(g, , (1) 
where Vi^s is the residual (particle- hole) interaction and Ilo{q,co) is the free polarization propagator: 

/rPk r Hi,)' 
(2^ J ^Go(|k + q|,u; + u;'))Go(k,u;') . (2) 

Go{k,co) is the free single-particle Green' s function: 

Goik,.)= '^'-'-^ + '^'^-'^ (3) 
' ' E-Ek + irj E-Ek-iv 

where Ek is the single-particle energy, kp the Fermi momentum and 9{x) the Heaviside step function, which equals 
1 for a; > and zero otherwise. The explicit expression for the free polarization in the nonrelativistic limit can 
be found in textbooks (see for instance Ref. [15]) and is usually referred to as the Lindhard function. Taking the 
low-momentum limit of this function, we obtain the Vlasov (semiclassical) approximation, whose accuracy increases 
for more relativistic particles [16]. In the present work, the nonrelativistic approximation is adopted for the nucleons, 
while the ultrarelativistic and semiclassical (Vlasov) approximation is taken for the polarization propagator of the 
electrons. 

The RPA equations describing the motion of a multifluid system take the matrix form: 

U"'{q,Oj) = Ul,{q,Oj) {5ik + vlijq)n^'^(q,uj)) , (4) 

which is an X iV matrix if TV particle species are present. The branches of the dispersion relation correspond to 
the poles of the response function. Disregarding the imaginary part of the free polarization, they are obtained by 
searching the zeros of the determinant A defined by: 

A = det [Re (1 - nov.es)] (5) 

where IIo and Vj-cs arc the frcc-polarization and residual-interaction matrices. It is a general feature of the RPA 
equations that each fluid generates a pair of branches if the interaction between its constituents is repulsive enough, 
and the interfluid coupling is not too strong. Typically, the pair of branches associated with the particle species i is 
positioned on each side of the line lu = qVpi , where Vpi is the Fermi velocity. This line corresponds to the singularity 
of the semiclassical Lindhard function (proportional to the free response function IIq in the Vlasov approach): 



L{si) = 2 - Sjln 
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with Si — Lu/{qVFi)- The lower branch of each pair is strongly damped. To liigliliglit this structure and help identify 
the nature of the different branches, the straight lines u = qVpi are drawn on all the multifluid dispersion relations we 
present in this work (dotted lines). To determine the relevant modes of excitations, it is essential to take into account 
the imaginary part of the response function: the quantity — Imll" gives the strength function associated with the 
particle species i, which shows peaks when collective modes arc present. For each excitation mode, the damping rate 
is measured by the width of the corresponding peak; the mode is actually a resonance if this width is not too large. 

III. PURE COULOMB INTERACTION 



Before considering the whole coupling scheme among neutrons, protons and electrons in neutron star matter condi- 
tions, we found it quite instructive to consider the proton-electron Coulomb coupling alone, disregarding the nuclear 
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interaction. This idealized case, which of course does not correspond to actual physical conditions, will permit clarifi- 
cation of two points, also relevant for the general physical situation. On one hand, it will make more transparent the 
role of electron screening on the proton dynamics and, on the other hand, it will give clear indications on the strengths 
of the electron and proton motion at the different excitation energies. In this way we will also describe in some detail 
the method we are going to use throughout the article for analyzing the excitation spectrum of homogeneous neutron 
star matter. 

As already mentioned in the Introduction, we are considering normal nuclear matter, leaving the case of superfluid 
matter to a future work. The range of baryon density will be fixed between the saturation density po = 0.16 fm^'^ 
and Spo: where indeed neutron star matter is expected to be homogeneous and not affected by possible "exotic" 
components like hyperons. Of course the proton fraction has to be fixed to choose the proton density. This is 
determined by /3 equilibrium and depends on the nuclear matter equation of state (EOS). This will be taken from 
the microscopic calculations of Ref. [17, 18], based on the many-body theory within the Bethe-Brueckner-Goldstone 
(BBG) scheme. The calculated P equilibrium takes into account the muon onset: however, for simplicity, muons and 
electrons will be treated in this work as a single ultrarelativistic fluid. We will work at zero temperature, but the 
results can be easily extended to finite temperature. The proton fraction as a function of baryon density is reported 
on Fig. 1. It is also given in Tab. I for selected values of the baryonic density, together with the corresponding Fermi 
momenta and electron plasma frequency (at pe — pp). Throughout the article, we use the unit system h = c = 1. 
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FIG. 1: (Color online) The proton fraction as a function of total baryon density in neutron star matter. The f3 equilibrium 
obtained from microscopic calculation is taken as a reference. The curves obtained with the Skyrme-type forces are shown for 
comparison. 
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0.16 3.7 327.2 110.5 6.15 

0.32 8.6 405.2 184.1 10.25 
0.48 13.0 456.2 242.3 13.48 



TABLE I: Density dependence of the proton fraction, Fermi momenta and electron plasmon frequencies in neutron-star matter 
conditions. The values given for the proton fraction Yp are obtained by microscopic calculation; these values are adopted for 
the calculations performed with all models presented in this work. 

For a one-component plasma, i.e. a charged Fermi gas on a rigid uniform background (jelly model), the RPA has 
been studied in detail, both in the nonrelativistic [15] and relativistic [19] cases. In this case, the residual interaction 
is Ui-cs = ''^c = Aire^ /q'^. The excitation spectrum, i.e. the line in the energy- momentum plane where the real part 
of the polarization function (1) has poles, is known to have a "thumblike" shape. This is illustrated on Fig. 2 (full 
line), where the case of a proton gas at a density p = po is considered. The upper branch appearing on this figure is 
called the plasmon excitation, which has a simple classical interpretation [20]. It is characterized by a finite energy 
at vanishing momentum: loq = Vf\/ (47re^)A^o/3, where Vp is the Fermi velocity and Nq is the level density. The 
lower branch is actually strongly damped and does not correspond to a real excitation that is able to propagate. This 
point will be further elaborated in the sequel. Prom the figure it is clear that in any case there is a maximum value 
of the momentum above which the spectrum stops and no excitation is possible. This is mainly a quantal effect, 
which is not present in the semiclassical approximation of RPA (Vlasov approximation). This is illustrated on the 
same Fig. 2, where the RPA and the Vlasov approximation spectra are compared. Apart from the momentum cutoff, 
one can see that the spectra are quite similar. The results for the relativistic case are analogous [19], except that the 
"thumb" becomes thinner and thinner as the particles at Fermi level become more and more relativistic [16]. In the 
ultrarelativistic case the two branches become extremely close near the cutoff, and the difference between RPA and 
Vlasov spectra appears smaller and smaller, apart again for the presence of the momentum cutoff. 

Having summarized the well-established results for the excitation spectrum with the pure Coulomb interaction 
within the jelly model, let us consider the system of protons and electrons, treating the two components on equal 
footing (pe model). In this case, the RPA equations are simply a 2 x 2 system, where the coupling is provided by the 
attractive Coulomb interaction between protons and electrons. The polarization propagator n**^ is now a 2 x 2 matrix: 
in Eq. (4) the indexes k stand for electrons (e) and protons {p). It has to be stressed that no exchange is included, 
i.e. only density-density correlations are considered. Since electrons are expected to be faster than protons, they are 
able to follow their motion and are more effective in screening the proton-proton interactions. The total spectrum is 
now formed by four branches, but only one is expected to have the characteristic property of a finite excitation energy 
at vanishing momentum: this will be identified with the plasmon mode of the proton-electron coupled system. The 
original proton plasmon mode is now reduced to a soundlike mode, since, due to the electron screening, the effective 
proton-proton interaction becomes of finite range (screened Coulomb interaction). This effect is illustrated by the 
short-dashed line of Fig. 2. This is also a well known phenomenon in classical plasma [20]. The screening can be 
easily seen in the present case by solving Eq. (4) for the proton-proton polarization propagator. In the pe model, we 
consider only the Coulomb interaction, and have the relation v^^ = v^^ = —v^'^ = Vc- After a short algebra, one finds: 



Since the unperturbed electron plasmon mode is substantially higher than the corresponding proton one, in the small 
q limit one can expand the free electron polarization propagator at zero frequency and one gets: 



This is equivalent to the RPA equation for a single fluid with a screened (finite-range) Coulomb interaction. For 
q < qc, the collective mode can be only a sound mode with a typical linear dependence on momentum. The inverse 
A — 1/qc is the electron screening length. One finds that q^ = 3(a;oe/VFe)^) where woe is the electron gas plasmon 
frequency at zero momentum and Vpe the Fermi velocity. In neutron star matter conditions, qc <C kpp and the 
screening length is much larger than the average distance between protons. Therefore, at increasing momenta the 
proton excitation should merge in the proton plasmon mode with no screening. The actual calculations confirm these 
expectations, as illustrated on Fig. 3 for three proton densities corresponding to total baryon densities po, 2po and 
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FIG. 2: (Color online) The proton dispersion relation in the jelly model. The energy cooe is the electron plasmon energy at zero 
momentum. The full line is obtained in the RPA scheme, the long-dashed line corresponds to the Vlasov approximation. The 
short-dashed line shows the effect of electron screening on the (Vlasov) proton plasmon branch. 



Spo- The reported branches are obtained by searching, at given momentum, the energies for which the determinant 
A of matrix Re [1 — IloVres] vanishes. In the pe model, A takes a simple expression: 

A = 1 - v,Re [ng + ng] (8) 

One observes that three branches go to zero at vanishing momentum. The upper branch is the plasmon mode of the 
electron-proton system, while the "thumblike" shape of the original proton branch in the uncoupled system can be 
identified with the lower " loop" . 

It has to be stressed that all branches, except the plasmon one, are actually damped. The polarization propagators 
Ho that appear in Eq. (4) are indeed in general complex. The best way to estimate the damping rate of each branch is 
to calculate the strength function — Imn*'(q, lo) (which is also a basic quantity for calculating the rate of many physical 
processes). In the present case, we consider the proton and electron components (i — p,e). At a given momentum q, as 
a function of lj, the strength function is expected to display peaks associated with the different branches, their width 
giving a measure of the damping. The electron and proton strength functions are reported on Fig. 4 at saturation 
density for selected values of the momentum q. One can see that peaks are present at the energies corresponding 
to the branches of Fig. 3, even if some slight shifts occur. Such shifts are expected, due to the contribution of the 
imaginary part. The proton and electron strength functions display peaks approximately at the same energy. The 
different heights at the peak maximum give a measure of the relative relevance of each components. The plasmon, up 
to a certain momentum, has no width, and the strength function presents a delta function singularity at the plasmon 
frequency. To put in evidence the plasmon branch, we have reported just a sharp bar at the corresponding energy. 
The relatively narrow peak at lower energy corresponds to the original proton plasmon in the uncoupled system. As 
discussed above, in the coupled system, electron screening reduces this branch to a soundlike mode; the coupling with 
electrons also produces a substantial damping. For this mode, the proton component is larger, but electrons follow 
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FIG. 3: (Color online) The four branches of the electron-proton system interacting by coulomb force only, under neutron star 
matter conditions. The energy uoe is the electron plasmon energy at zero momentum. Only the thick branches are associated 
with collective modes, while the thin ones are strongly damped. The straight lines u — qVpi that structure the spectrum are 
indicated with dots. 
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the proton oscillation, which allows an important screening at low momentum. However, as the momentum increases 
the screening is less effective, in agreement with the disciission above. On Fig. 4, we can observe that with increasing 
momentum (cases g = 5 — 25 MeV) the height of the peak increases and its width decreases. This tendency ends with 
the branch associated with the mode: at g = 30 MeV, where the branch has stopped, the peak is still present but low 
and wide. Let us finally comment the two other branches present in the spectrum (namely the lower branch of each 
pair). In the strength functions, they are associated with bumps that can be hardly identified with collective modes: 
they rather represent uncorrelated particle-hole excitations with a large spread in strength. Indeed, once the peaks 
are subtracted, these bumps resemble the ones present in a free Fermi gas. In any case, they can correspond only to 
strongly damped modes. The behavior at higher density is quite similar. 

An interesting check of the calculation is the fulfillment of the energy-weighted sum rules, which in the present case 
can be written: 



where wgi is the unperturbed plasmon frequency at zero momentum, for electrons {i = e) and protons {i = p). Notice 
that wop/woe ~ kp/M, where kp is the proton and electron Fermi momentum and M the proton mass. The sum rule 
is exact in the non-relativistic limit and valid also in the relativistic case in the limit of low momenta and energy, i.e. 
q <^ kp, Lu <^ Ep, which corresponds to the Vlasov approximation. The sum rules (9) arc satisfied by the calculated 
strength functions with an accuracy better than 10~^, which gives a check of the numerical accuracy. The contribution 
of each peak to the corresponding sum rule can be used to measure its relevance for many physical processes. Let 
us consider in particular the plasmon mode of the coupled system. Its contribution comes from a delta singularity, 
and can be calculated by standard numerical methods. The fraction of the sum rule for the electron component 
{i = e) due to the plasmon excitation is reported on Fig. 5 as a function of momentum. We see that the plasmon 
contribution decreases with increasing momentum. One can also observe that the electron strength scales as q/Kp 
for q/Kp > 0.07. This can be easily understood by noticing that in this region the plasmon energy u> ^ q (see Fig. 3), 
and then the electron polarization function depends only on q/Kp, since no other scale is present. For the proton 
component (i =1?), one can find that the contribution of the plasmon to the corresponding sum rule is never greater 
that few percentages and vanishingly small at the higher momenta. This confirms that the plasmon mode of the pe 
system is mainly an electron excitation. On the contrary, the lower mode, corresponding to the sound excitation, is 
mainly a (damped) proton excitation. 

The results we present in this article are obtained by solving explicitly Eq. (4). However, we found it interesting 
to check the validity of the static-electron approximation of Eq. (7) , noticing that at high density the proton plasma 
frequency is a substantial fraction of the electron one. In the strict static limit the proton strength function has a 
sharp peak (Dirac delta) corresponding to the soundlike mode. We found that the position of the peak is close to 
the maximum of the proton strength function corresponding to the damped mode obtained in the full calculation. A 
better approximation can be obtained by treating only the real part of Ilg in the static limit and leaving the imaginary 
part evaluated at the correct momentum. In this case, the approximate proton strength function around the sound 
mode is close to the exact one, as illustrated on Fig. 6. This is equally true for the electron strength function, which is 
also shown on the figure. Of course, the static-electron approximation cannot be kept for the whole energy range and 
in general one has to consider the exact strength function: the full response fmiction given by the matrix (4) has to be 
calculated in its complete form. This is what we will do in the rest of the article and is in line with previous studies 
on neutrino interactions in neutron-star matter [2-10]. Our discussion on the static approximation was intended only 
to assess its validity and limitation. 



The realistic physical situation in neutron stars includes the effective nuclear particlc-holc (p-h) interaction among 
neutrons and protons. The effect of p-h correlations on neutrino mean free path has been analyzed in a few articles [5, 
21]. Here we want to study in detail the spectral functions and the related collective modes, including the dynamical 
electron screening. We compare the results from effective (Skyrme) and microscopic nuclear models. We restrict 
the analysis to density fluctuation modes, which are expected to include the most collective ones. The particle-hole 
(or residual) interaction is derived in the Landau monopolar approximation. It amounts to considering only the Fq 
Landau parameter, neglecting the parameters Fl for L > 0. This is equivalent to averaging the effective interaction 
over the Fermi surface and is a reasonable approximation, provided the nonlocality of the effective interaction is not 
too large. This procedure not only simplifies the calculations, but also allows a more direct comparison between 
the phenomenological treatment of the Skyrme type and the microscopic one, where reliable EOS are available up 
to few times saturation density. Furthermore, it makes a clear connection between the residual interaction and 
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IV. ELEMENTARY EXCITATIONS IN THE REALISTIC CASE. 
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(Bare) 

Micro 


(0) 
-46.2 


(0) 

311.5 


(0) 
-394.1 


(1) 

0.926 


(1) 

0.856 


(0.349) 

0.376 


(0.118) 

0.137 


(4.05) 

3.75 


(1.37) 

1.17 


NRAPR 


47.1 


232.5 


-263.4 


0.806 


0.610 


0.432 


0.193 


3.27 


0.83 




SLy230a 


-172.8 


544.7 


-117.9 


0.544 


0.969 


0.640 


0.121 


2.21 


1.33 




(Bare) 
Micro 


(0) 
168.4 


(0) 
284.3 


(0) 
-131.7 


(1) 
0.897 


(1) 
0.760 


(0.432) 
0.481 


(0.196) 
0.258 


(5.02) 
4.50 


(2.28) 
1.73 


NRAPR 


252.1 


381.1 


155.4 


0.657 


0.447 


0.656 


0.438 


3.30 


1.02 




SLy230a 


88.3 


606.1 


317.8 


0.386 


0.870 


1.118 


0.225 


1.94 


1.98 


0.48 


(Bare) 
Micro 


(0) 
309.9 


(0) 
270.5 


(0) 
-5.2 


(1) 
0.951 


(1) 
0.848 


(0.486) 
0.510 


(0.258) 
0.304 


(5.65) 
5.37 


(3.00) 
2.54 


NRAPR 


398.9 


482.0 


459.6 


0.543 


0.357 


0.894 


0.722 


3.07 


1.07 




SLy230a 


320.2 


669.9 


588.8 


0.306 


0.746 


1.587 


0.345 


1.73 


2.24 



TABLE II: Density and model dependence of the residual interaction elements (bulk=momentum-independent part), effective 
masses, Fermi velocities and level densities. The bare values correspond to the free Fermi gas. The proton fraction Yp is fixed 
by the /9 equilibrium; the Yy value obtained by microscopic calculation is taken for all models. 



the thermodynamical properties related to the (free) energy-density curvature. Indeed, in the Landau monopolar 
approximation, the residual interaction reads [22]: 

where W is the energy density and A/o, = m*fcFj/7r^ is the level density, involving the effective mass m|; the indexes 
stand for neutrons (n) and protons (p). An equivalent way to express it involves the single-particle energy levels: 

= + Ui{k,Pn,Pp) , (11) 

where single-particle potential Ui{k) is obtained by the functional derivative of the interaction energy density with 
respect to the single-particle density pi. In the Landau monopolar approximation, we can identify the residual 
interaction with the derivative of the single-particle potential at constant momentum, taken at Fermi level: 

KL= (j^ikFi,pn,Pp)j ■ (12) 

It should be stressed that the nuclear p-h interaction in neutron-star matter condition is not well known. Only 
qualitative results can be derived, which are relevant to understanding the possible structure of the spectra and the 
connected phenomena that occur in neutron stars. In the framework of effective models such as Skyrme-Hartree-Fock, 
the residual interaction can be analytically deduced from the second derivative of the energy-density functional. We 
have obtained it applying directly Eq. (10) to the Skyrmc functional (sec Rcf. [22]). However, this fimctional is still 
poorly constrained for the high isospin asymmetries and densities involved in compact stars. The adopted Skyrme 
parametrizations are SLy230a [23] and NRAPR [24]. Both are modern Skyrme forces, whose fitting procedures take 
into account microscopic calculations of the neutron-matter EOS. However, important differences will be observed 
between these two models. On the other hand, from the microscopic point of view, we do not have an analytical 
expression for energy density 'H{pn, Pp)- The equation of state can be fitted by smooth functions, but the consequences 
of the fits on the curvature are not well controlled, especially at high asymmetry. Furthermore, it is also important 
to note that the k dependence of the single-particle potential Ui, which is connected to the presence of an effective 
mass, implies that even the interaction part of the energy-density curvature S'^Hint/Spf diverges for low values of pi, 
as explained in Appendix A. As a consequence, the isospin dependence of Hint cannot be treated in the parabolic 
approximation for high values of the asymmetry. For these reasons, we have chosen to extract the microscopic 
residual interaction rather from Eq. (12), performing numerical derivatives of Ui{kpi, p„, pp). The Brueckner-Hartrec- 
Fock (BHF) potential is calculated following the scheme of Ref. [17, 18], using the Argonne vis two-body interaction 
with Urbana IX three-body force. The BHF calculations also provide the effective masses to be used in the RPA 
equations. These effective masses must be carefully taken into account in the derivation of the residual interaction 
(see Appendix A). In all calculations the proton fraction, i.e. the nuclear matter asymmetry, is taken from the 
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FIG. 4: (Color online) Electron and proton spectral functions at saturation density and for selected values of the momentum. 
The vertical dotted lines show the energies at which the determinant A of Eq. (8) vanishes. The plasmon mode coincides with 
the vertical line of highest energy. The energy ojoe is the electron plasmon energy at zero momentum. Top: whole spectrum 
at q=10 MeV. Bottom: evolution of the protonlike mode with the momentum. The diagonal elements of the spectral function 
matrix are shown: — Imll'''' (full line), and — Imll'^'^ (dashed line). 
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FIG. 5: Fraction of the electron sum rule concentrated in the plasmon mode. The three curves correspond to baryon density 
equal to saturation density (diamonds), twice saturation density (circles) and three times saturation density (squares). 



corresponding microscopic calculations in neutron stars. This means that the calculations with Skyrme forces are 
not fully self-consistent. However, we do that to have a fair comparison between the Skyrme and microscopic results, 
since wc want to single out the structure of the excitation spectra under the same physical conditions, which should 

be the most realistic ones. 

The residual interaction elements w*^^ obtained with the different models are shown on Fig. 7 as a function of the 
baryonic density. The subscript b means " bulk" : this is the part of the residual interaction that is independent of the 
transferred momentum q, contained in the nuclear part. Note that, in the case of Skyrme forces, the nuclear residual 
interaction c;ontains a contribution proportional to q^; however, its effect is negligible at the relatively low q values 
considered in this work. In the following, the subscript b will be omitted. The effective masses and Fermi velocities 
are shown as a function of baryonic density on Fig. 8. Furthermore, all the model-dependent quantities of relevance 
for the present study are given in Tab. II, for the three selected values of the density we will consider. We see that 
the Skyrme models SLy230a and NRAPR produce quite different effective p-h interactions. Concerning the effective 
masses m*, we can note the atypical behavior of SLy230a, predicting a very low neutron effective mass (in particular, 
in contrast with the NRAPR and microscopic case, it predicts m* < m*). This results in a fast breakdown of the 
nonrelativistic description of the neutron fluid at increasing density, as is obvious from the superluminal values of the 
Fermi velocity that are reached even before 2po. Note that the NRAPR Fermi velocities also take large values as we 
approach 3po- Conversely, the quite high effective masses obtained in the microscopic framework are more consistent 
with a nonrelativistic treatment of the nucleons. 

We now turn to the study of collective density excitations including the nuclear force. The particle-hole nuclear 
interaction results in the coupling of protons with neutrons, and the possible appearance of neutronlike branches; 
furthermore, the effective mass it produces modifies the density level and Fermi velocities, affecting respectively the 
extension and orientation of the dispersion-relation branches in the {q, oj) plane. 
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FIG. 6: (Color online) Efltect of the static-electron approximation on the protonlike mode. The energy ojoe is the electron 
plasmon energy at zero momentum. The diagonal elements of the spectral function matrix are shown: — Imll'''' (with triangular 
marker), and — Imll*'^ (with round marker). 



Let us first study the simplified case where only the neutron and proton dynamics are considered (np model). The 
electron dynamics is suppressed, i.e. the electrons provide only a static uniform background of negative charge. Of 
course this is not a reahstic physical situation, but it will help the interpretation of the results of the more complete 
and realistic calculations. The protons interact also by the bare Coulomb potential and can display a plasmon mode. 
However, the nuclear interaction, in particular the coupling between protons and neutrons, can produce a shift and 
a width on the proton plasmon mode, as well as other modes where the nuclear effective interaction is essential. The 
dispersion-relation branches of this neutron-f proton model are displayed on Fig. 9 for the different nuclear models, at 
density 2po- In addition to the protonlike "thumb", two neutronlike branches appear around the line uj ~ VpnQ in the 
NRAPR and microscopic case, for which v™ is strongly repulsive; this is not the case with SLy230a, for which v™ is 
weaker. Since the line lu = VpnQ crosses the plasmon mode, in presence of the neutronlike modes an exchange occurs 
around this crossing: the upper branch, which at low q corresponds to the plasmon mode, turns into a neutronlike 
mode, while the plasmon continues on the upper part of the thumb. Note that in all three cases the protonlike mode 
is separated in two parts: at low momentum, where the frequency of the plasmon branch is above uj — Vpnl, the 
plasmon mode is undamped (Dirac delta); after the plasmon branch is crossed by the line lu — VpnQ, the imaginary 
part of the neutron response function induces a damping of the plasmon mode, which becomes a peak of finite width. 
This is shown by the strength functions displayed on Fig. 10 for momenta g = 10 and 20 MeV. We can see that, even 
when the plasmon mode is damped, it corresponds to a neat resonance peak with a quite small width; furthermore, 
it remains essentially a proton mode. The behavior of the neutron strength function deserves some comments. At 
q = 10 MeV, it does not present any peak, even in the NRAPR and microscopic cases where the relation dispersion 
presents neutronlike branches. These branches do not correspond to physical modes in this part of the spectrum, due 
to a strong damping. However, at g = 20 MeV where the highest branch is neutronlike, the NRAPR and microscopic 
forces show a " genuine" undamped mode (Dirac delta) near the edge of the neutron continuum. This mode is due to a 
strongly repulsive neutron-neutron interaction. A crucial role is also played by the effective masses. For the SLy230a 
force the neutron effective mass is smaller than for the other two forces, while for the proton effective masses it is the 
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FIG. 7: (Color online) Elements of the residual-interaction matrix (bulk=momentum-independent part) as a function of the 
baryonic density in neutron star matter, for the different nuclear models: microscopic approach and Skyrme forces. In the 
three cases, the proton fraction at j3 equilibrium is fixed by the microscopic calculation. 



contrary. For this force the neutron spectral function is expanded to a wider energy range, while the proton one is 
shrunk to a smaller range. 

All these findings are quite useful to understand the overall structure of the excitation spectra in the realistic 
case when the electrons are dynamically included in the calculations. The branches obtained for the full neu- 
tron+proton+electron (npe) system are shown for the different models, for three values of the baryonic density: 
po (Fig- 11), 2po (Fig. 12), and 3po (Fig. 13). The proton plasmon is shifted down by screening to a soundlike 
mode, as was the case without neutrons. This shift in turn modifies the neutron-proton coupling, in general reducing 
it, in particular for the SLy230a force, just because of the reduction of neutron phase space. However, as we have 
seen, at higher momenta {q > 10 MeV) the screening disappears and the scheme is the same as without electrons, 
except of course for the genuine plasmon mode at higher frequency that is almost completely an electron excitation. 
Proton-neutron excitations appear approximately decoupled from electron ones. Concerning the nuclear interaction, 
comparing the different spectra we observe that two kinds of quantities are crucial in determining their features (as 
was already observed in the np model): the neutron- neutron residual interaction v^^, and the effective masses. Thus, 
the neutronlike branches appear when the neutron-neutron interaction is sufficiently attractive: for all models, they 
are absent at density po, and appear at higher density. On the other hand, the effective masses affect two important 
quantities: the Fermi velocities Vpi — kpi/m* and the level densities iVpi = kpim* / {tt^). Since each pair of branches 
associated with a given particle species tends to follow the corresponding line lj — qVpi, neighboring Fermi velocities 
can lead to branch merging. This is what happens with NRAPR at two and three times po because of the high 
Fermi velocities of neutrons and protons, in contrast with the microscopic case (Figs. 12 and 13). The effect of the 
Fermi velocity is still more drastic with SLy230a, where Vpn takes unphysical superluminal values: this leads to the 
crossing of the electron plasmon branch by the line lj = qVpn- In this case the plasmon is strongly damped through 
the Coulomb interaction with protons and in turn through the coupling of the protons with neutrons by the nuclear 
effective interaction. However, this situation is not realistic and only shows the breakdown of the nonrelativistic treat- 
ment of the neutrons, inconsistent with their low effective mass. The second effect of m* concerns the level density: 
indeed, the presence of a pair of branches associated with the particle species i is favored by high values of Noi. This 
effect can be observed in different ways. Let us consider, for instance, the NRAPR and microscopic spectra at density 
2po. Despite a more repulsive interaction between protons, the neutronlike branches extend to higher momenta due 
to the higher neutron level density. Furthermore, although is more repulsive for NRAPR, the protonlike branches 
have a greater extension in the microscopic case because of the higher m* value (leading to a higher Nqp). We can 
also compare the SLy230a and microscopic neutronlike branches at density Spo; although is more repulsive for 
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FIG. 8: (Color online) Efltective mass (left) and Fermi velocity (right) as a function of the baryonic density in neutron star 
matter, for the different nuclear models: microscopic approach and Skyrme forces. In the three cases, the proton fraction at (3 
equilibrium is fixed by the microscopic calculation. 



SLy230a, the corresponding branches are vanishing due to the low m'^ value, in contrast with the microscopic case. 

We now study the strength functions — Imll" {i — n,p,e) obtained with the different nuclear models, for selected 
values of the momentum (g = 5, 15, 30 MeV): they are shown on Figs. 14, 15, and 16, respectively, for the baryonic 
density po, 2po, and 3po- 

Let us first comment the behavior of the proton mode. It is quite similar to the case of the pe model with pure 
Coulomb coupling: the mode is always damped, and is mainly a proton mode. At low momentum, the electron strength 
function in the peak region is close to the proton one (electron screening), especially at high density; however, at 
higher momentum the proton component of the peak increases, and becomes much more important than the electron 
one. The proton mode is suppressed in two kinds of situations. When the momentum exceeds the cutoff value of the 
protonlike branches, the proton peak becomes a wide bump, which cannot be identified any more with a collective 
mode. This is observed at saturation density, for q = 30 MeV in the NRAPR and microscopic cases (see Fig. 14). 
The other case of suppression occurs for NRAPR at 3pQ (Fig. 16), due to the merging of the branches (Fig. 13). This 
effect is due to the low proton effective masses: this corresponds to a high Fermi velocity Vpp, which induces a strong 
coupling with the neutron and electron components. Finally, we can notice the case of NRAPR at 2pQ and q = 5 
MeV. The proton mode appears in the strength function, and is quite similar to the typical one (Fig. 15); however, 
it is not associated with any branch of the dispersion relation, since at this value of momentum the upper proton 
branch has not yet appeared (see Fig. 12). 

Let us now consider the neutron strength function. In the realistic cases, the imaginary part of ImllQ is high 
when the proton peak occurs, and although the neutron strength function is affected it does not take the shape of a 
collective mode. A real neutron peak associated with the proton one is only obtained with SLy230a at high density, 
when Vpn is nonphysically high. However, a specific neutron mode occurs in some cases, usually indicated by the 
upper neutronlike branch. When it appears, this neutron peak is very sharp an high (exceeding the scale chosen for 
the present figures). We should note, however, that the presence of the neutronlike branches are neither necessary 
nor sufficient for the presence of such peak. Let us consider, for instance, the case of density 3po (Figs. 13 and 16). In 
the case of NRAPR at g = 30 MeV, a sharp neutron peak is obtained although there is no corresponding branch. On 
the contrary, in the case of SLy230a at g = 5 MeV, the upper neutronlike branch is not associated with a collective 
mode; the same kind of behavior has been observed above with the np model. 
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FIG. 9: (Color online) Dispersion relation in neutron+proton matter, at twice the saturation density and with proton fraction 
under neutron star matter conditions, for the different nuclear models. The energy cjoe is the electron plasmon energy at zero 
momentum. Only the thick branches are associated with collective modes, while the thin ones are strongly damped. The 
straight lines co = qVpi that structure the spectrum are indicated with dots. 
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FIG. 10: (Color online) Spectral functions in neutron+proton matter for selected values of the momentum, at twice the 
saturation density and with proton fraction under neutron star matter conditions, for the different nuclear models. The energy 
ujoe is the electron plasmon energy at zero momentum. The diagonal elements of the spectral function matrix are shown: 
— Imll"" (thin full line), and — Imll'''' (thick full line). The vertical dotted lines show the energies at which the determinant 
A of Eq. (5) vanishes. 



V. CONCLUSION 



The present work is dedicated to the study of normal density modes in homogeneous npe matter at /? equilibrium. 
Such modes should affect the thermal properties of neutron stars. We have stressed the behavior of the proton mode: 
although it corresponds to a plasmon mode in the jelly model, the inclusion of electron dynamics reduces it to a 
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FIG. 11: (Color online) Dispersion relation in neutron+proton+electron matter, at saturation density and with proton fraction 
under neutron star matter conditions, for the different nuclear models. The energy uoe is the electron plasmon energy at zero 
momentum. Only the thick branches are associated with collective modes, while the thin ones are strongly damped. The 
straight lines u = qVpi that structure the spectrum are indicated with dots. 
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FIG. 12: (Color online) Dispersion relation in neutron+proton+electron matter, at twice the saturation density and with proton 
fraction under neutron star matter conditions, for the different nuclear models. The energy ljob is the electron plasmon energy 
at zero momentum. Only the thick branches are associated with collective modes, while the thin ones are strongly damped. 
The straight lines u = qVpi that structure the spectrum are indicated with dots. 
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FIG. 13: (Color online) Dispersion relation in neutron+proton+electron matter, at three times the saturation density and with 
proton fraction under neutron star matter conditions, for the different nuclear models. The energy woe is the electron plasmon 
energy at zero momentum. Only the thick branches are associated with collective modes, while the thin ones are strongly 
damped. The straight lines ui = qVn that structure the spectrum are indicated with dots. 




FIG. 14: (Color online) Spectral functions in neutron+proton+electron matter for selected values of the momentum, at satu- 
ration density and with proton fraction under neutron star matter conditions, for the different nuclear models. The energy ujoe 
is the electron plasmon energy at zero momentum. The diagonal elements of the spectral function matrix are shown: — Imll"" 
(thin full line), — Imll*'*' (thick full line), and — Imll'^'^ (dashed line). The vertical dotted lines show the energies at which the 
determinant A of Eq. (5) vanishes. 



soundlike mode, as was clearly shown in the analysis of the pe system with pure Coulomb couphng. This effect is also 
observed for the full npe system, where the nuclear interaction that couples protons and neutrons is included. 

The nuclear quantities of relevance for the study of the density modes in npe matter (particle-hole interaction 
and effective mass) were found to be extremely model dependent: even the two modern Skyrme forces lead to very 
different predictions. This is essentially because of the high isospin-asymmetry values imposed by the /3 equilibrium 
in our study: for such neutron-rich matter, the energy density functional is still poorly constrained. The effective 
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FIG. 15: (Color online) Spectral functions in neutron+proton+electron matter for selected values of the momentum, at twice 
the saturation density and with proton fraction under neutron star matter conditions, for the different nuclear models. The 
energy uioe is the electron plasmon energy at zero momentum. The diagonal elements of the spectral function matrix are shown: 
— Imll"" (thin full line), — Imll^^ (thick full line), and — ImH'^'^ (dashed line). The vertical dotted lines show the energies at 
which the determinant A of Eq. (5) vanishes. 



mass was found to have an important effect on the dispersion relation obtained in the (q, uj) plane. It determines (i) 
the Fermi velocities that control the orientation of the different branches, possibly provoking some merging, and (ii) 
the level densities, which are an important quantity for the q extension of the nucleonlike branches. Note also that for 
a given repulsive value of uJ™, the neutron level density is crucial for the existence or not of the neutronlike branches, 
favored by high values of A^'on. It should be noted that low values of the nucleon effective masses at high density are 
not consistent with a nonrelativistic treatment. This is a serious limitation of the Skyrme models: in particular, with 
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FIG. 16: (Color online) Spectral functions in neutron+proton+electron matter for selected values of the momentum, at three 
times the saturation density and with proton fraction under neutron star matter conditions, for the different nuclear models. 
The energy cjoe is the electron plasmon energy at zero momentum. The diagonal elements of the spectral function matrix 
are shown: — ImH"" (thin full line), — Imll'''' (thick full line), and — Imll'^'^ (dashed line). The vertical dotted lines show the 
energies at which the determinant A of Eq. (5) vanishes. 



SLy230a super luminal values of the neutron Fermi velocity are reached even before 2po. In any case, we can stress the 
need to better constrain the isospin evolution of the nuclear effective masses: this could be allowed by experimental 
studies of exotic nuclei to be performed in the future radioactive-ion-beam facilities, for instance by the study of the 
isospin evolution of the level density. On the other hand, theoretical progress could improve the determination of the 
nuclear residual interaction, possibly through the extraction of Landau parameters with microscopic methods. 

An essential conclusion of our study is that the dispersion relation branches in the (q, uj) plane are not sufficient to 
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identify collective excitations in neutron-star matter: they have only an indicative role. To determine the existence 
and estimate the damping of such modes, we have to look for the peaks present in the strength function related to 
the imaginary part of the polarization matrix. Other contributions to the imaginary part, which we have neglected, 
come from two-particle— two- hole excitations [25] . However, our aim was to check for each mode if it was overdamped, 
and to this purpose the RPA approximation gives the correct answer. Considering the diagonal elements — Imll" 
with i = n,p, e, we have analyzed the density and momentum evolution of the different modes and identified which 
components are involved. In addition to the undamped electronlike plasmon mode, two nuclear modes have been 
observed. At high enough density, the neutron-neutron interaction becomes repulsive enough to induce a neutron 
mode, which appears as a sharp peak in the neutron strength function. Due to the strong Coulomb repulsion 
between protons, a protonlike mode corresponding to the screened proton plasmon is almost always present. Its 
strength increases with q until the region of the proton cutoff momentum, beyond which the mode is suppressed. It is 
interesting to note that a high value of the proton Fermi velocity can lead to the suppression of this mode by merging 
into the neutron continuum (case observed for NRAPR at high density). The proton collective mode is expected to 
play a role in the neutrino-emission process by pair breaking and recombination (PBC) in superfluid matter: this 
effect will be discussed in a future work. 
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APPENDIX A: RESIDUAL INTERACTION AND ENERGY-DENSITY CURVATURE 

In the Landau monopolar approximation, the bulk residual interaction (i.e. the part that is independent of the 
transferred momentum of a density fluctuation) is related to the curvature of the density energy H with respect to 
the single-particle densities pi (i = n.p). In this relation, the role of the effective mass has to be carefully taken into 
account. In particular, wc show in this appendix that the parabolic approximation for the asymmetry dependence of 
the interaction energy breaks down for neutron-rich matter, due to the presence of an effective mass. 

The energy density is the sum of a purely kinetic term K, and an interaction term Hint- In the mean- field approach, 
it is a functional of the single-particle density, from which are derived the individual energy levels. The energy level 
of a particle of species i = {n,p), bare mass rrii and momentum k takes the general expression: 

^2 

<^i{k,Pn,Pp) = + Ui{k,Pn,Pp) , (Al) 

where Ui is the mean-field potential. In general, this potential depends on k, which leads to the definition of an 
effective mass. The first derivative of the energy density corresponds to the Fermi level: 

= eFi{Pn,Pp) (A2) 

^2 

= et{kFi,Pn,Pp) = + Ut{kFi,Pn,Pp) ■ (A3) 

The second derivative reads: 

Introducing the effective mass m*{pn,Pp) defined by: 

^ikFi,Pn,Pp)) =^ (A5) 



and identifying the residual interaction 

(n^ n^^ = ( 

.5pj 



vlis{pn,Pp) = ( Yr{^Fi,Pn,Pp) ) , (A6) 
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we obtain the relation: 



Separating the interaction part of the energy density, we have: 
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(A8) 



Thus, if the effective mass is different from the bare one, the residual interaction does not simply reduce to the density 
curvature of Hint- This is due to the k dependence of the mean-field potential. Indeed, Eq. (A5) can be explicited as: 
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(A9) 



which is the standard definition of the local efi^ective mass at Fermi level. 

Let us now consider the proton-proton residual interaction {i = j = p) in neutron-rich matter. At low proton 

density, the curvature S'^TC-mt / 5 Pp diverges as 1/kpp because of the effective mass. Treating the isospin dependence of 
Hint in the parabolic approximation amounts to ignoring this efi'ect, and leads to a large underestimation of v^^^ in 
the case m* < nip. 
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